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We present a microscopic first-principles study of the neutral structural domain walls (DWs) in the multi- 
ferroic hexagonal manganites, which have been shown to exhibit cross-couplings between ferroelectricity and 
structural antiphase. We find that, in contradiction with previously proposed models, the interlocked antiphase- 
ferroelectric domain walls have approximately zero width, and their energy is lower than that of antiphase-only 
or ferroelectric-only domain walls. Furthermore, we show that the ferroelectric-only and antiphase-only DWs 
are superpositions of interlocked antiphase ferroelectric (AP+FE) DWs and intermediate domains inevitably 
exist through the DWs. Our results shed light on the question of why only AP+FE DWs are observed and how 
the topological defects emerge in polar hexagonal manganites. 
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In ferroic materials such as ferroelectric s, regions that dif- 
fer in the orientation of the ferroic order parameter (in ferro- 
electrics, the electrical polarization) are called domains. The 
planar defects between adjacent domains are called domain 
walls (DWs). In general, DWs have a different symmetry 
and structure than the bulk of the domains, and play an im- 
portant role in the switching mechanism of the primary or- 
der parameter. In multiferroic materials, which contain mul- 
tiple simultaneous ferroic orders, the DWs are of particular 
interest. For example, the ferroelectric DWs in multiferroics 
have been shown to exhibit novel properties such as local con- 
duction EHD and ferromagnetism [5 ]. In addition, couplings 
among multiple (anti-)ferroic orders, that are not observed in 
the bulk systesm, have been reported at the DWs (5HE1- 

The hexagonal manganites h-7?Mn03 (R=Sc, Y, Dy-Lu), 
with their coexisting ferroelectricity (7c ~ 1200-1500 K) and 
antiferromagnetism (rN~70-130 K), are currently among the 
most intensely investigated multiferroics [6, 9-11]. They are 
improper geometric ferroelectric s, in which the primary order 
parameter is a structural tilt trimerization of the MnOs poly- 
hedra, which is driven by minimization of the electrostatic 
potential lfTTtil3l . The trimerizaton corresponds to the con- 
densation of a zone-boundary K 3 mode of the high symme- 
try P63 jmmc structure, and lowers the symmetry to the polar 
P63 cm space group 1 12 ] but does not itself introduce a net fer- 
roelectric polarization. The trimerization can be described by 
two angles: the angle between the z axis and the Mn-apical O 
bond gives the magnitude of the tilting and can be regarded 
as the order parameter of the K 3 mode, and an azimuthal an- 
gle, <p, which describes the orientation of the tilting. Below 
the phase transition (p adopts one of six values separated by 
60°, which correspond to trimerization around three possible 
origins, with in- or out- orientation of the tilting; as a result 
six structural domains emerge, as shown in Fig0Ja) |fT4l . The 
polar T~ mode emerges as a result of an anharmonic coupling 
with the K3 mode, and the orientation of the subsequent ferro- 
electric polarization is set by the in- or out- orientation of the 
K 3 tilting. Then, <p = 0°, 60°, 120°, 180°, 240°, and 300° are 
often labeled a + , J3~ , y + , oT , J3 + , and y~ (Fig.[TJa)), where a, 
j3 and y represent different origins for the trimerization, and + 



or - indicate the out- or in- orientation of the tilts and corrse- 
pondingly the up- or down- orientation of the ferroelectric po- 
larization. 

There are six domains in 7?Mn03, and in principle, three 
types of structural DWs might be expected - ferroelectric 
only (FE), antiphase only (AP), and antiphase plus ferro- 
electric (AP+FE) - however, recent experiments showed that 
only AP+FE DWs exist, and that FE only DWs and AP only 
DWs do not occur [03] 02) . A net ferromagnetic signal has 
been reported at the AP+FE walls [5], as well as a wall- 
mediated coupling between the antiferromagnetic and ferro- 
electric domains (6). Finally, the intersections of the six pos- 
sible AP+FE DWs form characteristic vortex defects which 
have been shown to be topologically protected (U El 021 El- 
US). 

In this work we present a detailed systematic study of DW 
structures and energetics in YMn03 using first-principles den- 
sity functional theory. Our results are consistent with and 
complementary to a recent microscopy study using scanning 
transmission electron microscopy (STEM) fT9ll , which is not 
able to image the positions of light ions. In particular, we 
are able to explain the origin of the peculiar AP+FE selection 
rule. 

We begin by introducing our nomencalutre for the neutral 
DWs in h-7?Mn03. We first consider DWs between a + and 
/3~ domains. We define the (210) DW to be the wall of which 
normal vector from an a + domain to a f$~ domain points the 
[210] direction with respect to the primitive unit cell vectors 
(Fig.[l{b), left wall). The (210) DW is then the wall with the 
same orientation, but the a + and J3~ domains are reversed, as 
shown on the right of Fig. [TJb). Although the detailed struc- 
tures are not known yet, the structures of the (210) and (210) 
DWs are in fact different, as can be seen for example by com- 
paring the distances between the layers of R ions staying at 
2a Wyckoff position |20j|i n each domain, £ 2a ^ as the wall is 
traversed. We see in Fig. [lib) that iia-^j-' (3w 3 +1) fl (where n 

is an integer) across a (210) DW and l%a — ^ ' ^J^ a across 
a (210) DW. Since the geometries are different the structures 
and energetics must be in principle different. On the other 



2 



(a) OY • lower O • upper O 
V® ^ ^ > © ^ > 

^-y^-© \ \ * © \ 

©-■*©- ®— (5 ©*■ s -{•) 

a+(<p=0°) a-((p=180°) 



(b) 



„ _V3 (3n+1) a , „ _V3 (3n+2) a^ 
"» 3 2 3^^ 




(210) DW ^ ®, ^.(210)DW * ® ^ * 

.©-«©-<•--& 

j8 ' rra ' 0+ 
2 a 



G> ^ © Q-^->-0 101 °J [210] 

«^g>«®-v v© -«*v [2io] [ 10 °] 
2a y ^ © * x v © >( ^ 

jS+ (a>=240°) jS- (a>=60°) © > X © © X y © © > X © 

/©V ^ / © V 

0..y....^...0 (g)-^-V-<g) <gj -*>-e*0(no)DW0 o ^ 0(110) dw£x) o -e>(g) 

* © O&afc © v x X ® » ® H X ® 

b ©% V ^ / * ^ ^ 

c y+ (<p=120°) y- (<P=300°) a+ 0" a+ 

„ _V3 (3n+2) a ^, - _V3 (3n+1) a, 
**~2 3 \ c 2«~2 3/*\ 

-> V ©> © 

((210)DW ^ © ^ H 

V-^® V ^ ^ © x 

~ (110) (8)-^-- : >-0) ,0- 



< d > (010) ( 1 20) (no) (e) 

(110) _ \ | / _(210) 

-►(100) 




(210) / 
(110) 



(120) 



(010) 




FIG. 1: (a) Schematics of the 6 possible structural domains in 
YMn0 3 . The view is along the z axis, perpendicular to the direc- 
tion of ferroelectric polarization. The big and small circles show re- 
spectively the positions of Y ions and adjacent oxygens at the apices 
of the trigonal bipyramids, and their arrows indicate the direction of 
their displacements from the high symmetry P6 3 /mmc phase. We la- 
bel each domain by its phase <p, defined to be the (counter)clockwise 
angle of the displacement directions of upper (lower) oxygens rel- 
ative to a + domain fT4l . Note that the 2a Wyckoff positions vary 
regarding the antiphase domains, (b-c) Our nomenclature for the 
neutral DWs sandwiched between a + and pr domains. Note that the 
(210) and (210) DWs are different whereas (110) and (110) DWs are 
same by symmetry (see text for details), (d) The classification of 
the equivalent DWs between a + and J3~ domains that are shown with 
normal vector directions, (e) DWs between a + and y~ domains. The 
(210) and (210) DWs are reversed from those between a + and /T 
domains. 



hand, the (120), (110) and (210) DWs (normal vector from 
a + to pr pointing [120], [110] and [210] directions) are all 
equivalent to each other by symmetry, and the (120), (110) 
and (210) DWs (defined analogously) also form a symmetry- 
equivalent set. We refer to these two types of DWs as {210} 
and {120} DWs respectively. The two DWs perpendicular to 
the [110] direction are (110) and (110) DWs (Fig. (TJc)), and 
they are identical by symmetry. (100), (010), (100), and (010) 
DWs are also equivalent to (1 10) and (110) DWs, and we refer 
to this type of DW as {110} DW. The DWs between p + and y~, 
and between y + and a~ are equivalent to that between a + and 
J3~, for the same DW direction. Furthermore, the DW between 
a~ and J3 + is equivalent to that between a + and J3~ as identified 
by m z symmetry operation. If f3~ is replaced by y~, however, 
the DWs are reversed as shown in Fig.[TJd). 

First principles spin-polarized calculations were done with 
projector-augmented wave method [2D as implemented in 
VASP l22l . The exchange-correlation interactions among 
electrons were treated by the local density approximation with 
Hubbard U correction |]23l[24]|. The parameters were set to 
the values of previous reports (U = 8 and / = 0.88 eV on the 



Mn-3d orbitals) l25lL and A-type magnetic configurations are 
adopted. Lattice constants for the calculation of DWs were 
fixed to the values of the relaxed unit cell without DWs, and 
internal positions were optimized in each case until the forces 
acting on all atoms converged to less than 0.005 eV/A respec- 
tively. We carefully tested the convergence of the plane wave 
cutoff energy and &-point sampling. The cell size dependence 
was checked by adopting several supercells containing up to 
300 atoms. 

We begin by calculating the energies of the two DW models 
originally proposed by Choi et al. in Ref. ifTSTl . These have the 
orientation of the {120} and {210} DWs in our nomenclature, 
and are obtained by averaging the structures of two neighbor- 
ing domains so that the Y ions along the wall retain their high 
symmetry positions, and the DWs retain two-fold rotation and 
oglide reflection symmetries. Therefore we refer to them as 
the 2/c models. These 2/c models are physically reasonable 
since walls that are the structural average of two different do- 
mains, and which therefore retain aspects of the paraelectric 
structure, are known to be stable in perovskite-structure fer- 
roelectrics [26-28 ]. Figure [2ja) shows our calculated low- 
est energy walls within the constraints of the 2/c symmetry 
for {120} and {210} DWs in a 300-atom supercell. First, we 
note that, after structural relaxation, only the constrained Y 
atom directly at the boundary adopts its paraelectric position; 
on either side the Y ions immediately adopt their bulk po- 
sitions, indicating that the DW is only one unit cell (u.c.) 
wide. Second, we note that the intermediate structure at the 
DW has (p « 30° . This is reminiscent of the centrosymmetric 
P3c InMn03 structure in which | of the In ions retain six- 
coordinated high-symmetry positions and the MnOs bipyra- 
mids trimerize at angles intermediate to those of the YMn03 
structure (<p = 30° + n • 60°) lfT3l . We find also that the tilting 
magnitudes of the MnOs bipyramids 6 are lower at the DWs 
since the tilting correlate with the displacement of Y ions lfT2l 
and | of Y ions keep staying at paraelectric positions at the 
DWs. Indeed, the tilting magnitudes of YMn03 in P6^cm and 
P3c structures are respectively 5.1 and 4.5 degrees. The 2/c 
DW energy, calculated with four supercells then linearly ex- 
trapolated to infinite spacing yields 390 meV/supercell =195 
meV/DW = 195 meV/(a • c) « 44.8 mJ/m 2 . This is close to the 
energy difference per unit cell between the P6^cm and P3c 
structures (155 meV/u.c.) |[T3lL consistent with the fact that 
the width of each DW is approximately 1 u.c. of P3c struc- 
ture. The slightly higher energy would come from the walls 
between the P6^cm and P3c structures. 

Next, we displace the centrosymmetric Y ions at the DWs 
along the z axis as shown in Fig. [2jb) and fully relax the in- 
ternal positions, without the constraint of the symmetry. Our 
subsequent structural optimization does not recover the 2/c 
DW models and instead yields lower energy DWs in which 
two-fold rotation symmetry breaks and two-fold screw rota- 
tion symmetry appears. In these low energy walls, <p and 
change abruptly across the DW from one bulk value to the 
next, and we can consider the width of {120} and {210} DWs in 
YMn03 to be effectively zero; we refer to them subsequently 
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FIG. 2: Calculated {120} and {210} DW structures for (a) the 2/c DW models and (b) the sharp models (see text). The average Mn-apical O 
distance, and tilting magnitude (6) and azimuthal angle ((f) of Mn0 5 bipyramids for each layer are also shown. Dashed lines indicate the bulk 
quantities, which are recovered at the middle of each domain. Two-fold rotation and c-glide reflection symmetries exist in the 2/c DW models 
whereas two-fold screw rotation and c-glide reflection symmetries exist in the sharp models. The uniform Mn-O ap distances even through the 
DWs indicates Mn0 5 bipyramids rigidly tilt even through the DWs. 



as sharp DWs. We note that our calculated {210} DW is con- 
sistent with recent STEM measurement of the Tm positions 
in hexagonal TmMn03 for walls within the bulk of the sam- 
ples (T9). It is worth pointing out that {210} DW for the 2/c 
model is observed at the edge of the samples in TmMn03. 
This may be due to the truncation of the electrostatic potential 
or defect accumulation at the edge of the sample. The wall 
width is much thinner than expected for a ferroelectric DW, 
which is typically a few u.c. wide, but rather characteristic 
for an antiphase boundary. The DW energy does not strongly 
depend on the supercell size (Fig. [5]), indicating that the ad- 
jacent walls are neither strongly repulsive nor attractive. Our 
calculated DW energy extrapolated to infinite spacing is 11.2 
mJ/m 2 , which is a quarter of the value calculated for the 2/c 
DWs. Microscopically, the DWs intersect the shared edges of 
the MnOs bipyramids and are staggered when viewed from 
c-axis (Fig. |2[b)). The length is 2/ V3 = 1.15 times longer 
than that of the straight line observed in mesoscopic exper- 
iments such as piezoresponce force microscopy (PFM) (see 
also Fig.[9ja)). In 7?Mn03 polar T 2 mode, i.e. ferroelectric po- 
larization couples to the K3 mode and the polarization is pro- 
portional to the amplitude of K 3 mode Q K3 at Q K3 > 0.5 (T2J. 
Since 6 is almost constant even close to the DWs and thus 
Qk 3 ~ 1, the polarization direction should abruptly reverse 
across the DWs with similar bulk polarization. 

The supercells considered above each contain both {120} 



and {210} DWs, and the calculated DW energy is the aver- 
age value. To extract individual DW energies, we next con- 
struct supercells containing only one type of DW. Since the 
structural phase (p shifts +60° and -60° in the [210] direc- 
tion through {210} and {120} DWs, respectively, a supercell 
containing six identical DWs is required to retain the periodic 
boundary condition. We use 60- and 240-atom supercells for 
{210} DW, and 120- and 300-atom supercells for {120} DW; 
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FIG. 3: Energies for {110}, {120} and {210} DWs calculated with dif- 
ferent supercell sizes. The dash-dotted lines are linear extrapolations 
to infinite spacing (zero inverse distance). Plus signs (+) indicate the 
averages of {120} and {210} DW energies calculated with the super- 
cells containing both DWs. As expected, these values correspond 
approximately to the average of the individual {120} and {210} DW 
energies calculated with the large supercells containing one type of 
DW each. 
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FIG. 4: Calculated supercells for evaludating the {210) and {120) AP+FE DW energies as well as structural phases (p for each layer, (a) 60- 
and (c) 240-atom supercells contain six {210) DWs and (b) 120- and (d) 300-atom supercells contain six {210} DWs so as to retain the periodic 
boundary condition. In all cases the DWs are evenly spaced. 
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FIG. 5: (a) Unit cell of YMn0 3 . Mn ions staying at different layer 
and Y ions displacing in the opposite direction are distinguished by 
different colors, (b) Calculated {110} DW structure for the sharp 
model with a 300-atom supercell. Only half of the supercell is shown 
with consideration for the supercell symmetries. Numbers indicate 
the average structural phases if at each column of Mn0 5 bipyramids. 
The thick line indicates the {110} DW derived from the phase shift. 
Note that the {110} DW is regarded as a combination of {120} and 
{210} DWs. 



in all cases the DWs are evenly spaced (Fig. [4]). Figure [3jb) 
shows our calculated {210} and {120} DW energies from the 
two different supercells. Surprisingly, although the distances 
between DWs in 60- and 120-atom supercells are only -\a 

and ~y -\a, respectively, the estimated DW energies are very 
close to those with 240- and 300-atom supercells. This in- 
dicates that interactions between the DWs are very weak. We 
find that the energy of the {120} DW is slightly higher than that 
of {210} DW. The Y ions displace up-down-up-down across 
{210} DWs and up-up-down-down across {120} DWs (Fig. [2}, 
which likely contributes to differences in the electrostatic en- 
ergies. 

Next we calculate the structure and energetics of {110} 
DWs. Figure [5] shows our relaxed {110} DW structure cal- 



culated within a 300-atom supercell. Based on our results for 
{120} and {210} DWs, we constructed the initial structure by 
displacing the Y ions along the z axis so that domains sud- 
denly change through the DWs. This sharp configuration is 
retained on fully relaxing the supercell, and our calculated 
DW width is effectively zero. We note that two-fold screw 
rotation and oglide reflection symmetries exist in the super- 
cells for this wall orientation, and so the {110} and {110} DWs 
are symmetry-equivalent and we can extract their energies di- 
rectly from the supercells containing both walls. Our calcu- 
lated infinite-separation DW energy extrapolated from calcu- 
lations for four supercells is 12.5 mJ/m 2 (Fig. [3}, which is 
slightly higher than those of {120} and {210} DWs In fact the 
{110} DW can be regarded as a combination of {120} and {210} 
DWs (Fig. [5]), and indeed the {110} DW energy estimated from 
the {120} and {210} DW energies is 2/ V3 • 11.2 = 12.9 mJ/m 2 
which is close to our explicitly calculated value. We note 
also that the length of the {110} DW viewed from c axis is 
4/3 = 1.33 times longer than that of the straight line. 

Our finding that the DW width is almost zero is in contrast 
to the behavior in conventional perovskite counterparts such 
as PbTi03, where the widths of 180° DWs are narrow but fi- 
nite (> 1 u.c. width) [26-28 ], and the ions in the wall region 
transition from one orientation to the other through their high- 
symmetry paraelectric positions. The adoption of the high- 
symmetry position can be regarded as a geometric frustration 
of the positions in the two opposite domain structures. This 
frustration does not exist in h-7?Mn03, because the up-up- 
down or down-down-up displacements of the R ions, deter- 
mined by the tiltings of the MnOs trigonal bipyramids, can be 
locally preserved through the DWs, as seen in Figs [2] and [5] 
Instead, the translational periodicity that is broken at the DWs 
is reminscent of stacking faults commonly observed in close- 
packed structures. It is also interesting to note that the 180° 
DW energy in YMn03 is considerably lower than those of 
perovskites with similar Curie temperatures, for example the 
corresponding 180° DW energy in PbTi0 3 (T c =l 65 K) is 132 
mJ/m 2 (26) and in BiFe0 3 (T c =l 100 K) it is 829 mJ/m 2 l27l . 
While in conventional ferroelectrics, Curie temperatures often 
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FIG. 6: (a) Schematics of winding AP+FE DWs, which are combi- 
nations of {120), {210), and {110} DWs. The dashed lines indicate 
the mesoscopically observed DWs and the directions of Y ion dis- 
placements are changed through the DWs. The thick lines indicate 
the DWs derived from the phase shift, (b) DW energy as a function 
of the DW direction, which is geometrically interpolated with {120}, 
{210}, and {110} DW energies represented as £{120}. £{2io}> and £"{110}, 
respectively. 



correlate with the magnitude of the ferroelectic polarization, 
this is not the case in YMn03 where the Curie temperature is 
high (1258+14 K) [20] and corresponds to the tilting transi- 
tion (I2l, while the improper ferroelectric polarization is not 
large (5.6 yuC/cm 2 ). 

Experimentally, most DWs do not form straight lines cor- 
responding to a single orientation, but wind, as illustrated in 
Fig.[6]for the case of our calculated lowest energy sharp walls. 
It is clear from Fig.[6](a) that the winding DWs consist of com- 
binations of {120} (green), {210} (red), and {110} (blue) DWs, 
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FIG. 7: Calculated (a) AP and (b) FE DWs perpendicular to the [210] 
axis with 270- and 300-atom supercells, respectively. The models 
locally preserve the up-up-down or down-down-up displacements of 
the R ions across the DWs, resulting in the interstitial domains. 



which in turn are combinations of {120} and {210} DWs. We 
can then estimate the energy of a DW pointing in arbitrary 
direction from geometrical relations. Our results are shown 
in Fig. [6jb), where we see that the {210} DW has the low- 
est energy. This is consistent with the recent observation of 
stripe domain patterns in h-7?Mn03 grown below the trimer- 
ization temperatures lfT8l . The stripes develop along the [110] 
direction, which could correspond to either {120} and/or {210} 
DWs. From our results, we suggest that the DWs are likely 
to be of {210} type. Since {210} DWs do not vanish when 
they meet together for the same direction, they should lead to 
a topological protection of the stripe domains under electric- 
field application even without the topological defects. 

We now perform structural relaxations for FE-only and AP- 
only DWs, to address the question of why such walls do 
not occur in h-7?Mn03. Based on the results of the AP+FE 
DWs, we constructed initial models to preserve the up-up- 
down or down-down-up displacements of the R ions through 
the DWs. As shown in Fig. [7] the relaxed structures retain the 
off-centering of the Y ions. By analyzing the tilt patterns of 
the MnOs polyhedra, and the corresponding Y off-centering 
patterns, we see that an AP DW in fact is formed from two 
AP+FE DWs, and a FE DW is composed of three AP+FE 
DWs. Indeed, the average DW energies that we obtain for AP 
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FIG. 8: (a) Initial structure for calculating the migration of {210} 
AP+FE DW with NEB method, (b) The calculated energy landscape 
through the DW migration. The dashed line indcates the average DW 
energy of {210} and {120} 2/c DWs. (c-e) The initial, intermediate, 
and final structures in the frame in (a). 



and FE DWs with 270- and 300-atom supercells are 23.0 and 
34.4 mJ/m 2 respectively; almost exactly two and three times 
the average value for the AP+FE DWs (11.2 mJ/m 2 ). This is 
because the phase difference between two domains with dif- 
ferent origin (polarization) is 120° (180°), whereas AP+FE 
DWs can change the phase only +60°. The obtained DW en- 
ergies are not consistent with the six- state clock model used 
in Ref. lfT8l for Monte Carlo simulation, in which the ratio of 
AP+FE, AP, and FE DW energies is assumed to be 1:3:4. We 
also see that the minimum possible width of the intermediate 
phase is ^ • \a or ^ • |a, depending on the type of the AP+FE 
DW. Jungk et al. have found a width of 60+10 nm for the in- 
terstitial domain with PFM after applying electric field |[T6ll . 
This would be due to the depolarization field and/or the reso- 
lution of PFM. Finally we note that, in our energy minimiza- 
tions, the AP and FE DWs do not tend to yield larger regions 
of the intermediate phase calculations; this is because, as we 
discussed earlier, the interaction between DWs is very weak 
and so there is not a strong driving force for them to sepa- 
rate. At the high temperatures close to the phase transition 
temperature, however, we expect AP+FE DWs to randomly 
distribute for entropic reasons, and the artificially narrow in- 
termediate regions shown in Fig. [7] would naturally broaden 
so that AP-only and FE-only DWs would not survive. 

In ferroelectric s an electric polarization is reversed by the 
application of an electric field through the migration of the 
DW, and thus it is valuable to understand the microscopic mi- 
gration mechanism. In Fig. [8] we present the energy landscape 
through the DW migration and the transition structure calcu- 
lated with the nudged elastic band (NEB) method l29l . It is in- 
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FIG. 9: (a) Experimentally observed domain pattern with PFM re- 
produced from Ref. [13] and (b-c) pictorial view of three Y sites and 
their possible displacements around the topological defect. Note that 
around the center, a +{ ~\p +{ ~\ and y +(_) share the same site, and thus 
two of them displace up (down) and one of them down (up). 



teresting that the transition structure and its energy at the mid- 
dle point are almost same with those of the 2/c DW (Fig. [2]). 
In other word, the 2/c DW is the local maximum of the mi- 
gration energy curve. This is because when the DW moves, Y 
ions specified in Fig.[8jc) go through the paraelectric sites, and 
simultaneously the MnOs trigonal bipyramids rotate. We also 
see that the DW width increases from approximately zero to 1 
u.c. in the migration, and finally moves from one edge of the 
intermediate 2/c DW to another, meaning the minimum mi- 
gration step is 1 u.c. The energy barrier height is 33.2 mJ/m 2 , 
and is almost three times of the DW energy. Although the bar- 
rier height is comparable to PbTiC>3 (37 mJ/m 2 ), the ratio to 
the DW energy is much higher than that of PbTiC>3, of which 
the energy barrier is only 36% of the DW energy l26l . How- 
ever, the energy barrier height would significantly reduce at 
the high temperatures close to the phase transition tempera- 
ture, and AP+FE DWs would result in random distribution. 

Based on our finding that the lowest energy AP+FE DWs 
have the sharp configuration, in this last section we discuss 
the likely structure of the topological defects that form at the 
meeting points between six AP+FE DWs. Figure [9] shows the 
displacements of Y ions in the AP+FE DWs near the topolog- 
ical defect. Of the six Y ions around the central point, half of 
them displace up and half down along the z axis. This causes 
a frustration of the Y ion at the center of the topological de- 
fect which can not displace in a direction that allows all DWs 
to maintain their sharp configuration. It is possible, therefore, 
that the Y ion at the center might remain at its paraelectric 
position; otherwise a modification of the wall structure would 
be required in the vicinity of the topological defect. We note 
that this picture does not correspond to that of Ref. fT5lL as the 
latter paper used the higher energy P2/c models for the DWs 
meeting at the defect. It is clear that further theoretical and 
microscopy studies would be desirable for a full understand- 
ing of the structure of the topological defects. 
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